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ABSTRACT 

We outline a method for fitting binary-lens caustic-crossing microlensing events based 
on the alternative model parameterisation proposed and detailed in Cassan (2008). As 
an illustration of our methodology, we present an analysis of OGLE-2007-BLG-472, 
a double-peaked Galactic microlensing event with a source crossing the whole caustic 
structure in less than three days. In order to identify all possible models we conduct an 
extensive search of the parameter space, followed by a refinement of the parameters 
with a Markov Chain-Monte Carlo algorithm. We find a number of low-% 2 regions 
in the parameter space, which lead to several distinct competitive best models. We 
examine the parameters for each of them, and estimate their physical properties. We 
find that our fitting strategy locates several minima that are difficult to find with 
other modelling strategics and is therefore a more appropriate method to fit this type 
of events. 

Key words: gravitational microlensing - data modelling - cxtrasolar planets - binary 
stars - robotic telescopes 
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1 INTRODUCTION 



[1931), 



Gravitational microlensing (jPaczvriskil Il986l ) occurs when 
the light from a source star is deflected by a massive com- 
pact object between the source and the observer, leading 
to an apparent brightening of the source, typically last- 
ing a few days to a few weeks. When the deflecting body 
has multiple components, such as a planet orbiting its 
host star, there can be perturbations to the brightening 
pattern of observed sources. These perturbations can be 
large even when caused by low-mass objects, making them 
detectable using small ground-based telescopes. Modelling 
these lightcurve anomalies can lead to the detection of sub- 
tle effects, allowing for measurements of prop erties such as 
the so urce star limb-darkening coefficients fe.g. lCassan et all 
20041). the mass o f stars with no visible companions (e.g. 



Ghosh et al] 120041). and the detect i on of extrasolar planets, 
as sugge sted bvllVIao fc Paczvriskil (|l99ll ) and first achieved 



in 2003 (|Bond et alj|2004 ). 

Nevertheless, anomalous microlensing events usually re- 
quire very detailed analysis for a full characterisation of their 
nature to be possible, making them challenging. This applies 
in particular to a class of microlensing events which display 
caustic crossing features in their lightcurves. These events 
are of primary interest, because they account for around 
ten percent of the overall number of detected microlenses, 
and they represent an important source of information o n 
physical properties of binary stars (jjaroszvnski et al"]|2006l ). 
However there exist several degeneracies that affect the mod- 
elling of this type of events: without a robust modelling 
scheme and a full exploration of the parameter space, it 
is impossible to pin down the true nature of a given event. 
In addition to this, calculations of anomalous microlensing 
models for extended sources are very demanding computa- 
tionally. 

Given these issues, brute force is not an option when 
modelling caustic-crossing events, and one has to devise 
ways of speeding up calculations, for example by exclud- 
ing regions of parameter space which cannot reproduce fea- 
tures that appear in data sets. A way to achieve this is to 
use a non-standard parameterisation of the binary-lens mod- 
els that ties t hem directly to data features, as proposed by 
ICassanl (|2008l ) . which we recall below. 

In this paper, we present our method for exploring the 
parameter space, and describe our approach to find all pos- 
sible models for a given event (Sec. [2}. We then use OGLE- 
2007-BLG-472, a microlensing event observed in 2007 by 
the OGLE and PLANET collaborations, as an illustration 
of our methodology applied on a binary lens event which 
intrinsically harbors many ambiguities (Sec. [3]). We finally 
discuss the implications of the individual competitive mod- 
els that we find in order to discriminate between realistic 
microlensing scenarii. 



2 BINARY-LENS EVENTS FITTING SCHEME 

2.1 Parameterisation of binary lens lightcurves 

A static binary lens is usually described by the mass ratio 
q < 1 of the two lens components and by their sep aration d, 
expressed in units of the angular Einstein radius (|Einsteinl 
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where M is the mass of the lens, and Dl and Ds are the 
distances to the lens and the source respectively. Such a lens 
produces caustics where the magnification of the source di- 
verges to infinity for a perfect point source. The positions, 
sizes and shapes of the caustics depend on d and q. For the 
binary lens case, caustics can exist in three different topolo- 
gies, usually referred as close, intermediate and wide; bifur- 
cation values between t hese topologies are anal ytical expres- 
sions relating d with q (jErdl fc Schneider!! 19931 ) . In the close 
regime, there are three caustics: a central caustic near the 
primary lens component, and two secondary caustics which 
lie off the axis passing through both lens components. In the 
intermediate case, there is only one large caustic on the axis. 
In the wide case, there is a central as well as a secondary 
caustic, both on the axis. The limits between these config- 
urations are i ndicated as th e dashed lines in e.g Fig. [2] (see 
also Fig. 1 of lCassanl (|2008l )). 

The description of the lightcurve itself requires four 
more geometrical parameters in addition to d and q. In the 
current standard parameterisation of binary lens lightcurves, 
these are the source trajectory's angle a with the axis of 
symmetry of the lens, the time of closest source-lens ap- 
proach to the binary lens centre-of-mass to, the Einstein 
radius crossing time £e and the source-lens separation at 
closest approach uq (in units of #e). Finally for a uniformly 
bright finite size source star, we add a further parameter, 
th e source size p * in units of #e- However, and as discussed 
in ICassanl (|2008l ). this parameterisation is not well adapted 
to conducting a full search of the parameter space, because 
the value of the parameters cannot be directly related to 
features present in the lightcurve, namely caustic crossings 
for the type of events we are discussing in this paper. Conse- 
quently, most of the probed models in a given fitting process 
do not exhibit the most obvious features in the lightcurve, 
leading to very inefficient mod elling. 

To avoid this drawback, ICassanl (|2008l ) introduced a 
new parameterisation in place of a,to,uo an d te which is 
closely related to the appearance of caustic crossing fea- 
tures in the lightcurve. The caustic entry is then defined 
by a date t Gn try when the source center crosses the caustic^ 
and its corresponding (two-dimensional) coordinate gentry on 
the source plane. Howeve r, since by de finition this point is 
located on a caustic line, ICassanl l|2008h introduced a (one- 
dimensional) curvilinear abscissa s which locates the cross- 
ing point directly on the caustic, so that (entry = C( s cntry). A 
given caustic structure is fully parameterised by < s < 2. 
The caustic entry is then characterised by a pair of param- 
eters (ientry, Sentry), and in the same way the caustic exit by 
(icxit, s ox it). These four parameters (in addition to d, q and 
p») which describe the caustic crossings therefore also define 
an alternative parameterisation of the binary lens, far better 
fitted to describing the problem at hand. 



Alternatively, any other point at a fixed position from the 
source center can be defined as a reference. 
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2.2 Exploration of the parameter space 

We start by exploring a wide region of the parameter space 
with a (d, q) grid regularly sampled in logarithmic scale. 
This choice comes from the fact that the size of the caustic 
structures behave like power-laws of the lens separation and 
mass ratio, and so do the corresponding lightcurve anoma- 
lies. We fit for the remaining model parameters tcntry, ^exit, 
Sentry, s cx i t and p, , (d, q) being held fixed. From this, we then 
build a x 2 (d, q) map that we use to locate the best-fit (d, q)- 
regions. As mentioned previously, there exist binary lens 
configurations which involve central and secondary caustics. 
In these cases (i.e. in t he wide and close binary cases) and 
following ICassar] l|200Sf ). we study separately models where 
the source crosses the central or the secondary caustic by 
building two X 2 (d,q) maps, corresponding to each configu- 
ration. 

In order to sample efficiently and extensively Sentry and 
Scxit (which determi ne the source trajec tory), we use a ge- 
netic algorithm (e.g. ICharbonnead[l995t ) that always keeps 
the best model from one generation to the next one (elitism). 
In fact, since we consider only models displaying caustics at 
the right positions, there are a couple of local minima asso- 
ciated with different (s en try, s ex it) pairs which would usually 
be missed by other minimisation methods, while a genetic 
algorithm naturally solves this problem in an efficient way. 
However, since such an algorithm never converges exactly 
to the best model, we finally refine the model by perform- 
ing a Markov-Chain Monte-Carlo (hereafter MCMC) fit: we 
start several chains and use the criterion bv lGewekd (|l992l ) 
to assess convergence to a stationary posterior distribution 
of the parameter probability densities. 

From the obtained \ 2 maps, we then identify all the 
local minima regions and use the corresponding best models 
found on the (d, q) grid as starting points to refine the pa- 
rameters, including (d, q) that we now allow to vary. Since 
the fit is performed within a minimum \ 2 region, the fitting 
process is very stable and fast. 



3 APPLICATION TO OGLE-2007-BLG-472 

3.1 Alert and photometric follow-up 

On 19 Augus t 2007, the OGLE Early Warning System 
(Udalski 2003) flagged microlensing candidate event OGLE- 
2007-BLG-472 at right ascension a 2 ooo.o = 17:57:04.34, and 
declination fcooo.o = -28:22:02.1 or I = 1.77°, b = -1.87°. 

The OGLE lightcurve has an instrumental baseline 
magnitude 1=16.00, which may differ from the calibrated 
magnitude by as much as 0.5 magnitudes. Lensing by the 
star in the point source-point lens (hereafter PSPL) approx- 
imation accounts for a broad rise and fall in the lightcurve, 
peaking around MHJE0=4334 .0 with an apparent half- 
width at half-peak of about 10 days (Fig. [1]). Although the 
observed OGLE flux rises only by 0.06 mag in the non- 
anomalous part of the lightcurve, the shape of the curve 
hints that blending is important for this target, with only ~ 
12% of the baseline flux due to the un-magnified source. 

On 19 August (MHJD=4331.5) an OGLE data point 
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Canopus 1-Orr 

Danish 1,54m 
OGLE 
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Figure 1. OGLE, UTas and Danish data set for OGLE-2007- 
BLG-472 data sets. Data points are plotted with 1-cr error bars. 
. The x-axis is time in HJD-2450000. 



Telescope 


Data 


Error bar rescaling factor 


UTas 1.0m 


34 


1.79 


Danish 1.54m 


84 


1.55 


OGLE 


857 


1.21 



Table 1. Datasets and error bar rescaling factors. 



showed sudden brightening of the source, with subsequent 
PLANET (UTas Mt. Canopus 1.0m telescope in Tasmania 
and Danish 1.54m telescope at La Silla, Chile) and OGLE 
data indicating what appears to be a fold caustic crossing by 
the source, ending with a PLANET UTas data point on Au- 
gust 21 (MHJD=4334.1). The caustic entry is observed by 
a single OGLE point, while the caustic exit is well covered 
by our UTas data set (Fig. [TJ . Treating the lightcurve as the 
addition of an anomaly to a PSPL lightcurve, the underly- 
ing PSPL curve then apparently reaches peak magnification 
on August 22 (MHJD=4335.45). Particularly crucial in our 
data set is the UTas observation taken within a few hours 
of the caustic exit, which puts constraints on the position 
of the caustic exit on the lightcurve, and on the size of the 
source. Although some V-band observations were taken, the 
V lightcurve of this event is flat and does not allow us to 
place constraints on the V flux parameters. 

3.2 Data reduction 

We reduced the PLANET data for this event using the data 
reduction pipeline pysis3.0 (Albrow 2008). This pip eline uses 
a kern el as a discrete pixel array, as proposed by iBramichl 
|200&t ). rather than a linear combination of basis func- 
tions. This has the advantage that it removes the need for 
the user to select basis functions manually, which can lead to 
problems if inappropriate functions are chosen. In addition 
to this, the pixel array kernel copes better with images that 
are not optimally aligned. The result of using this pipeline is 
a better reduction than was obtained with other methods. 
We kept all points with seeing <3.5 arcseconds. Although 
some dubious points remain with this simple cut, the size 
of their associated error bars reflects their lack of certainty 
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and ensures their weight in modelling procedures is appro- 
priately reduced. Our final data set consists of 34 UTas data 
points, 84 points from the Danish 1.54m telescope, and 857 
points from OGLE (Table [TJ. 

3.3 Modelling OGLE-2007-BLG-472 

After a first exploration of the parameter space, we find a 
best model (close to model C c , see below) which we use 
as a basis to rescale our error bars. In fact, these can vary 
rather widely from one telescope to another and are often 
underestimated by photometry software. Ignoring this effect 
would misrepresent the relative importance of the data sets. 
From this step, we choose the rescaling factors shown in 
Table[T] obtained by setting x 2 /d.o.f. ~ 1 for each data set. 
We then use the rescaled data to perform a new parameter 
space exploration. 

We then apply the fitting scheme detailed in section [5] 
to our data sets. In particular, we choose a spacing between 
the (d,q) grid points of 0.070 in logd and 0.275 in logg. 
For the genetic algorithm fit, we use a model population 
of 200 individuals evolving during 40 generations, which has 
proven to be enough to safely locate the regions of minimum 
X 2 ■ Finite source effe cts are comput ed using the adaptive 
contouring method of iDomini kl ([20071). 

The final X 2 (d,q) maps that we obtain are plotted in 
Fig. [2] for the intermediate and central caustic configura- 
tions, and Fig.[3]for the intermediate and secondary caustic. 
The red crosses show the underlying (d, q) grid, and the blue 
shaded contours indicate values of A\ 2 ~ 5, 20, 50, 100, 250, 
where the reference model is C„, the global best- fitted model 
(as obtained in Section [33} . 

3.4 Excluding minima 

Fig. 2] shows a zoom of this region (secondary caustic and 
close configuration d < 1), with an overplot of ts isocontours 
(orange lines) roughly equally spaced on a logarithmic scale. 
With this fitting approach, we put no initial constraints on 
the Einstein time £e, though it will always remain physical 
(£e > 0). Since we are not using any Bayesian prior for 
this parameter, we find that very good fits to the data are 
obtained with values of £e > 300 days, which correspond to 
the minimum region in the left lower part of Fig. [3] Such long 
Einstein times are not likely to happen commonly, and it 
may happen that some of the values found for to correspond 
to a lightcurve that reaches its peak well in the future; these 
are very unlikely to be acceptable solutions. Hence, instead 
of using a prior for £e i n the fitt i ng pr ocess, we adopt the 
posterior distribution of iDominikl (|2006h , from which we see 
that £e > 400 days, well in the tail of the distribution, can be 
used as a cut for a model to be physically plausible. Thus in 
the following, we will not consider solutions with values of £e 
greater than 400 days. This means that we will not include 
the low-g (q ~ 0.001) minima in the following discussion. 

Although a very well-covered lightcurve generally en- 
ables a good characterisation of the deviation caused by 
the caustic approach or crossing, degeneracies make find- 
ing a unique best - fittin g mo del difficult. In particular, 
iGriest fc Safizadehl (|l998h and IDominikl (|l999l ) identified a 
two-fold degeneracy in the projected lens components sep- 
aration parameter d, under the change d <-> 1/d, when 



q <C 1. Moreover, iKubas et al] (|2005l ) showed that very simi- 
lar lightcurves could arise for a source crossing the secondary 
caustic of a wide binary system and for the central caustic of 
a close binary system. These degeneracies cause widely sep- 
arated x 2 minima in the parameter space, which must then 
be located by exploring the parameter space thoroughly. In 
addition to these degeneracies, imperfect sampling can in- 
crease the number of local x 2 minima; short event in partic- 
ular are prone to under-sampling, leading to difficulties in 
modelling. OGLE-2007-BLG-472 is no exception, as shown 
in the next section. 

3.5 Refining local minima 

We see from Fig. [5] (intermediate and central caustic) that 
there are three broad local minima in the region around the 
white filled circles marked as C c , I and W c ("I", "C" and 
"W" for intermediate, close and wide models respectively, 
and subscript "c" for central caustic). In Fig. [3] (intermediate 
and secondary caustic), a best-fit region can easily be located 
around the region marked C s (subscript "s" for secondary 
caustic), besides region I. 

Now allowing for the parameters d and q to vary as well, 
we use our MCMC algorithm to find the best solutions in 
each of these local minimum regions. These are identified 
with white filled circles on Fig. [2] and [3] and correspond to 
the models listed in Table and shown in Fig. [SJ El El an d 
[8] The best model lightcurve is dominated by strong caus- 
tics, which all viable models must reproduce, with the low- 
magnification base PSPL curve barely noticeable. All mod- 
els have the first anomalous OGLE points on the descending 
side of the caustic entry except for the worst model, model 
W c , which has this OGLE point on the ascending part of the 
caustic entry. Statistically, the former case is more likely to 
be observed since the ascending part of the caustic entry 
happens much more rapidly than the descending side. 

Our best model, C s , has x 2 = 949 for 975 data points, 
with the other competitive models at A^ 2 = 13.2 (model 
C c ), Ax 2 = 23.5 (model I) and A X 2 = 39.6 (model W a ). 

3.6 Discussion 

Fig. 2] shows that the models with a source crossing a sec- 
ondary caustic have increasingly large values of £e as they 
go towards lower values of the mass ratio. This is expected 
since the time At between Gentry and i ex it is fixed by the data. 
As the size of caustics scales with q 1 ^ 2 , and ts ~ At/q 1 ^ 2 , 
the source must therefore cross the Einstein Ring over a 
longer timescale in order to conserve the right timing for 
Sentry and s xit- In addition to this, blending decreases for 
decreasing values of q , and therefore decreases with increas- 
ing £e, contrary to what might be expected. Indeed, one 
would expect the blending factor g — Fb/Fs (where Fb and 
Fs are the blend and source flux respectively) to increase 
with increasing tv in order to mask long timescales and re- 
produce the observed timescale. However in this region of 
parameter space, the caustics are weak, which means that 
too much blending would not allow models to reproduce the 
observed rise in the source magnitude at the caustic entry 
and caustic exit. For a region of parameter space to contain 
satisfactory models, there must therefore be a fine balance 
between blending, timescale and mass ratio. 
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Figure 2. x 2 (^>9) map for the intermediate and central caustic configurations. Contour lines and minima regions (in blue shades) are 
plotted at Ax 2 = 5, 20, 50, 100, 250. The two dashed curves are the separation between the close, intermediate and wide regimes. The 
models are labelled and marked with white filled circles. 
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Figure 3. Same as Fig. [2] for the the intermediate and secondary caustic configuration. 



For models where the source crosses a central caus- 
tic, the impact parameter uo must decrease with decreasing 
mass ratio, since the size of central caustic decreases with 
decreasing mass ratio, and the range of allowed itn decreases 
if the source must cross the caustic. This means that for 



smaller mass ratios, blending will have to increase in order 
to mask the correspondingly higher PSPL magnification of 
the source that results from the smaller impact parameter. 
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Secondary caustic, close separations zoom 




Separation (d) 



Figure 4. Map of the value of in the (d, q) plane for converged models at each grid point, superimposed on the x 2 map, zoomed in 
on the close regime part of parameter space. Contours lines (orange) are labeled with their corresponding value of ie while \ 2 contour 
lines are plotted at A\ 2 = 5, 20, 50, 100, 250 and filled with gradual shades of blue. The dashed curve is the separation between the 
close and intermediate regimes. The models of Table [2] are labelled and marked with white filled circles. 



3.7 Physical properties of the models 

3. 7. 1 Source characteristics 

A colour-magnitude diagram of the field (Fig. [9} was pro- 
duced extracting 1497 stars from I and V images at t = 
4340.08 (I) and t = 4340.13 (V) taken at the Danish 
1.54m telescope. The combination of the source and the 
blend lies very slightly blueward of the red giant clump, at 
(V — 7)=2.43. All the models, however, are heavily blended 
(Table [2]|. The actual source magnitude and blending mag- 
nitude for each model can be found using the equations 
la = /base + 2.5 log(l + g) and I b = I s — 2.5 log(p). 

Using this equation, we find source magnitudes ranging 
from 17.89 (model C B ) to 20.21 (model C c ) (see Table [2). 
Our V-band data set does not allow us to determine the 
source's colour, but assuming that the source is a main se- 
quence star we use the calculated I magnitude of the source 
for each model to est imate a colour, using the results of 
iHoltzman et alj (|l998t ) . This then enables us to estimate the 
source's angular radius which we use in Section ^. 7. 2l to com- 
pute probability densities of the lensing system's properties. 

We calibrate the baseline magnitude of our target 
(source and blend combined) using the location of the red 
clump as a reference. We find /base = 15.61 ± 0.10, which is 
in agreement with the OGLE value of /base = 16.00 ± 0.50. 
Comparing this to the location of the red clump, we can de- 
rive an estima te for the reddening c o efficie nt A\. From Hip- 
parcos results, IStanek fc Garnavichl (|l99&f ) find an absolute 
magnitude for the red clump at Mi,ho = — 0.23±0.03. Using 
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Figure 9. Colour - Magnitude diagram of the field. The target 
OGLE-2007-BLG-472 is shown as a black triangle at {V -1,1) = 
(2.43, 15.61). The position of the dcblendcd source for each model 
is labeled and indicated by a dotted line and a coloured diamond, 
with the blend for each model also plotted as a diamond in the 
same colour. 



a distance modulus to the galactic centre of n = 14.4 1 ±0.0 9 
(i.e. assuming Ds = 7.6 kpc) ([Eisenhauer et all 120051 ). 
this translates to a dereddened magnitude for this target 
of /base = 14.18 ± 0.09. Hence using the relation Ai — 
/base — A/i,ac — Mi we get a value for the /-band reddening 
parameter of Aj = 1.43 ±0.13. Alternatively, fitting 2MASS 
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Figure 5. Best-fitting binary lens model C s with residuals and a zoom on the anomaly (left inset). Data points are plotted with l-cr 
error bars. The trajectory of the source in the lens plane with the caustics is plotted as an inset in the top right corner of the figure, 
with the primary lens component located at the coordinate system's origin. 
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Figure 7. Same as Fig. [5] for model C' c 

Model W„ 



o 



O 
CM 



fi - 



Canopus 1.0m 1 
Danish 1.54m 
OGLE 




. s ' 




4315 4320 4325 4330 



4335 4340 



4345 



4350 



HJD-2450000 

Figure 8. Same as Fig. [5] for model W c . 
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Table 2. Best-fitting binary lens model parameters. The blending factor g(I) = Fb{I)/ F$(I) is given for the OGLE data (/-band). 
The error bars were rescaled for each telescope by the factor given in Table [T] which lead to the rescaled \ 2 indicated here. Physical 
parameters are also given for each model, for the case of a lens in the disk, and a lens in the bulge. These were calculated using the 
procedure detailed in Sec. 13.7.21 



Parameter 


ivioacl u s 


ivioacl o c 


ivioael l 


Ivioacl W c 


Units 


X 2 (rescaled a) 


949.00 


963.16 


972.48 


988.55 




Ax 2 


_ 


13.2 


23.5 


39.8 


_ 


*UTas 


23.79 


24.83 


26.41 


28.86 




-^Danish 


79.77 


79.60 


80.75 


88.93 




^OGLE 


845.50 


858.77 


865.24 


870.55 




to 


4587.18 ±0.80 


4332.27 ±0.29 


4332.10 ±0.27 


4334.99 ± 0.28 


MHJD 


'e 


213.82 ± 1.04 


52.00 ±3.63 


38.32 ±2.60 


53.46 ±0.81 


days 


a 


2.810 ±0.006 


3.227 ±0.030 


3.305 ± 0.037 


4.570 ±0.018 


rati 


UQ 


-1.573 ±0.013 


0.091 ±0.005 


0.164 ±0.019 


0.277 ±0.010 




p*/io- 3 


0.34 ± 0.01 


0.98 ±0.19 


1.55 ±0.16 


1.33 ± 0.05 




d 


0.427 ±0.002 


0.673 ±0.011 


0.760 ± 0.015 


2.158 ± 0.0169 




q 


0.078 ± 0.001 


0.177 ±0.017 


0.236 ± 0.024 


0.288 ± 0.0096 




g(I) = F B (7)/Fg(7) 


7.15 ±0.013 


68.11 ±0.013 


40.13 ±0.09 


56.98 ± 0.019 




h 


17.89 ± 0.01 


20.21 ± 0.01 


19.65 ±0.09 


20.02 ±0.01 




h 


15.75 ± 0.01 


15.63 ± 0.01 


15.64 ±0.09 


15.63 ±0.01 




(V-I)s 


1.80 ±0.10 


1.93 ±0.11 


1.91 ±0.11 


1.92 ± 0.12 




0* 


1.18 ±0.24 


0.46 ±0.09 


0.59 ±0.12 


0.50 ±0.10 
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isochrones to our CMD, we obtain a value Ai = 1.46 ± 0.08 
and E(V — I) = 1.46 ± 0.11. We use these values of red- 
dening to determine dereddened magnitudes and colours for 
the source of each model. These, together wi t h the surface 
brightness relations from iKervella fc Fouquel (|2008h , allow 
us to calculate the apparent angular radius of the source 0, 
for each of the models, given in Table [2] 



3.7.2 Lens characteristics 

Although the characteristics of any microlensing event de- 
pend on various properties of the lensing system, including 
the mass of the lenses, the only measurable quantity that 
can be directly related to physical properties of the lens is 
the timescale of the event tE. While the physical properties 
of the lensing system can be fully constrained when the pho- 
tometry is affected by both finite source-size effects and par- 



allax, when these are not measured, such as is the case with 
our analysis OGLE-2007-BLG-472, we can still use Bayesian 
inference to determine probability densities of physical prop- 
erties of the lens, based on a chosen Galactic model. We have 
chosen not to include parallax in our analysis because its ef- 
fect would be very small for such a low-magnification event; 
in addition to this, we are only seeking a first-order analy- 
sis of binary-lens events with our current method, although 
second-order effects such as parallax and lens rotation will 
be taken into account in future work. 

We use our fitted value of the source size parameter p* 
to place constraints on the mass of the lens, which can be 
expressed as a function of fra ctional distanc e x = Dl/Ds 
and the source size p* as (e.g. |Pomiriik1ll99cf ) 

M{x) _ c 2 D s Ql x 
Mq 4GM q pi 1 - x ' U 
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Figure 10. Probability densities for the mass of the primary lens star and the fractional distance Dj^/Dq, for a lens in the disk (left 
side) and a lens in the bulge (right side). The values quoted in Tables [2] & [2] are the median value and the limits of the 68.3% confidence 
interval. On each plot, the probability densities are plotted for model C s (red), model C c (green), model / (dark blue), and model W c 
(light blue). 



where M is the mass of the lens, f?* is the angular radius 
of the source, the value of which is given in Table [2] and 
other quantities are denned as before. The mass-distance 
curve showing constraints from this equation is plotted on 
Fig. ED 

However, since we cannot measure parallax for this 
event, w e use a probabilistic approach following that of 
IDominikl (|2006l ) to derive probability densities for physical 
properties of lens components. Th e Galactic mode l used here 
is a piecewise mass spectrum (e.g. IChabrier||2003i ) . two dou- 
ble exponentials for the disk mass density and a barred bulge 
tilted a t an angle of 20° with the direction to the Galactic 
centre (|Dwek et alj|l995h . and the distribu tion of effective 
transverse velocities used in IDominikl (|2006l ). 

Using these galactic models, we infer properties for the 
lensing system, separating the cases where the lens is in 
the Galactic disk and in the Galactic bulge. For a lens in 
the disk, we find a primary mass 1.50~^q ^Mq and a sec- 
ondary mass of O.12lg'QgM0, at a distance of l.OOto'-jj; kpc 
with a lens velocity of 25lg 4 kms -1 . For a lens in the bulge, 
we find a primary mass 41 j^Mo and a secondary mass of 



3.2^\'\Mq, at a distance of 6.7^q ' 7 kpc with a lens veloc- 
ity 1671^° kms -1 . These are the physical lens properties for 
the lowest-x 2 model (model C s ). The values of these phys- 
ical parameters for the other models are given in Table [2] 
Probabilty densities of these properties for all models are 
plotted on Fig. [TO] 

3.7.3 Discussion 

For our lowest-^ 2 model, the parameters we find imply very 
unusual properties of the lensing system. As discussed in 
Sec. 13.41 the fact that we find these types of models is a con- 
sequence of the fitting approach we are taking. Traditional 
fitting methods would struggle to find these minima, since 
most of them require providing a starting point in parameter 
space. This is an issue when solely using an MCMC algo- 
rithm to fit microlensing events: although an MCMC run 
may be able to make its way through parameter space to 
find minima reasonably far away from its starting point, it 
is highly unlikely that a chain will be able to reach a mini- 
mum that has parameters different from the starting point 
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Figure 11. Mass-distance diagram showing the constraint on the 
lens mass from the source size, given by Eq. J2), for each model. 
The curves are labeled with the name of the model to which they 
correspond. 



by more than one order of magnitude. As we see from Fig. [4] 
there exist minima in many parts of parameter space, with 
values of £e that are different by almost two orders of magni- 
tude. These parameters are non-intuitive, since they cannot 
be guessed only by looking at the lightcurve. As a result, it 
is improbable that this kind of parameters will be used as 
starting points for "classic" fitting algorithms. 

We solve this problem for the static binary-lens case 
by resorting to the method described in Sec. 12.21 Using this 
approach, we manage to systematically locate minima in 
parameter space. However, one then has to be careful with 
interpreting the significance of the obtained model param- 
eters. The shape of probability densities shown in Fig. [10] 
for model C s indicates that our value of £e push the lens 
mass towards the end of the adopted mass spectrum in the 
Galactic model we have adopted. This results in the abrupt 
transients seen on Fig. [10] Similarly, the mass-distance curve 
for model C s on Fig. [TT] shows that the mass of the lens in- 
creases very rapidly for lenses above ~ 1 kpc. These unusual 
curves are caused by a value of £e ~ 200 days. Models with 
£e ~ 3000 days (corresponding to the law-q minimum visible 
on Fig.[3]&:[4]| are obviously not acceptable, but how can we 
formally reject them? Finding these models from minima in 
the x 2 surface shows the limits of using \ 2 as a strong crite- 
rion for favouring models. A solution to this would be to use 
prior distributions on as many of the parameters as we can. 
During the MCMC part of our fitting process, this would 
mean that we obtain posterior distributions that are differ- 
ent from the ones obtained without using prior distributions 
on the parameters, or, equivalently, assuming uniform pri- 
ors for all parameters. Such priors can be obtained in var- 
ious ways, such as looking at the distribution of timescales 
for past microlensing events or calculating th ese distribu- 
tions from Galactic models (e.g. iDominild [20061 ) . or by using 
luminosity functions of the G alactic bulge to find a prior 
for the blending factor g (e.g. iHoltzman et al.lll998l ). Such 
work requires careful consideration of which priors are most 
appropriate to use, and is beyond the scope of this paper. 
Using these priors in combination with our method to find 
minima will lead to more robust determination of minima 



by taking into account our knowledge of physical parameter 
distributions. 



4 SUMMARY AND PROSPECTS 

Our analysis of OGLE-2007-BLG-472 is a good illustration 
of the importance and power of using parameters that are re- 
lated to actual observed features. Indeed, despite incomplete 
coverage of the caustic entry and high blending, a few cru- 
cial data points and an appropriate choice of non-standard 
parameters enable us to find several good binary-lens model 
fits to our data for this event by exploring the parameter 
space systematically. Some of the good fits that we identify 
have unphysical parameters, and we must then reject them. 
However using this parameterisation allows us to be certain 
that the parameter space has been thoroughly explored. We 
find four models with different parameters: two close binary 
models, one intermediate configuration, and a wide binary 
model. The lowest-;^ 2 model corresponds to a G dwarf star 
being lensed by a binary system with component masses 
1.50tJ;||M© (for the primary) and 0.12tg;^M© (for the 
secondary), which are compatible with our blending values. 
However it is obvious from physical parameter distributions 
that using \ 2 as a s °l e criterion for determining the best 
model is insufficient, because it does not take into account 
our knowledge of the distributions of physical parameters. 

Since the approach presented in this paper can form 
the basis for a systematic, wide ranging exploration of the 
parameter space to localise all possible models for a given 
data set, it is particularly relevant to current efforts to au- 
tomatise real-time fitting of binary-lens events. This could 
prove useful to provide faster feedback on the events be- 
ing observed and prioritise observing schedules, especially 
on robotic telescopes. Expanding robotic telescope networks 
controlled by automated intelligent algorithms are expected 
to play an increasingly impor tant role in microlen sing sur- 
veys in the coming years (e.g. iTsapras et al]|200ct ). Fitting 
methods such as the one described in this paper are essen- 
tial for making sure any anomalies are interpreted correctly, 
and that minima are located in as large a part of parameter 
space as possible. 
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